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\& ■ Abstract 
O 

t^- ' In this work we discuss two methods of calculation of quark propagation in the quark-gluon 

O ' 

plasma. Both methods make use of the Nambu-Jona-Lasinio model. The essential difference of 

o : 

these calculations is the treatment of deconfinement. A model of confinement is not included in 



the work of Gastineau, Blanquier and Aichelin hep-ph/0404207 , however, the meson states they 
consider are still bound for temperatures greater than the deconfinement temperature T c . On the 
other hand, our model deals with unconfined quarks and includes a description of the qq resonances 
found in lattice QCD studies that make use of the maximum entropy method (MEM). We compare 
the qq cross sections calculated in these models. 
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It is found in studies of heavy-ion interactions at RHIC that the system very rapidly 
reaches equilibrium and that a model based upon hydrodynamics is appropriate for the early 
stages of the collision in which one expects to form a quark-gluon plasma [1]. Recently 
Gastineau, Blanquier and Aichelin have proposed a model of critical opacity based upon 
their calculation of the qq interaction for temperatures in the range of to 350 MeV [2]. 
A central feature of their model is the temperature dependence of masses of the u, d and 
s quarks, as well as the masses of the tt and K mesons. (See Fig. 1 of Ref. [2].) In their 
model of the qq scattering amplitude they calculate s, t and w-channel exchange of it and K 
mesons in the temperature range < T < 350 MeV. (We remark that deconfmement takes 
place at T c ~ 170 MeV in unquenched QCD calculations and at T c = 270 MeV in quenched 
calculations.) In Ref. [2], the matrix elements in the s and t channels are given by 

-iM t = S CUC3 5 c ^ C4 u(p 3 )Tu(p 1 )[i'Df(p 1 - Pz)]v{p4)Tv(p 2 ) (1) 
+ S CuC3 6 C2tCi u{p 3 )(i'y 5 T)u{p 1 )[iVf{p 1 - Pz)\v{p±){i-f b T)v(p 2 ) 



-iM s = S cl>C2 S C3;C4 v(p 2 )Tu(p 1 )[iV^(p 1 +p 2 )]v(p4)Tu(p 3 ) (2) 
+ 5 CuC2 5 c ^ c JJ(p 2 )(i^ 5 T)u(p 1 )[iV^(p 1 + p 2 )]v(p 4 )(ii 5 T)u(p 3 ), 

where Pi(p 2 ) is the momentum of the incoming q(q) and ^3(^4) that of the outgoing q(q). 
The Cj are color indices and the various T's are the isospin projectors on the meson states. 
Here. V s and V p are the meson propagators of the form obtained in the NJL model, 

and 

l-G P J p) {) 

with J s and J p being the polarization tensors in the scalar and pseudoscalar channels, 
respectively. (In Eqs. (3) and (4) we have modified the notation of Ref. [2] to be in closer 
correspondance to the notation we have used in our work.) 

In Ref. [2] the following Lagrangian is used without the term proportional to Gy 
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C = q(ij - m°)q + -f ^[(gA^) 2 + (qi^Xq) 2 } (5) 

i=0 

" ^E[(^7^) 2 + (^757^) 2 ] 

+ ^{det[g(l + A 5 )g] + det[g(l - X 5 )q}}. 

Here, m° is a current quark mass matrix, m° = diag(m^,m^,m^). The Aj are the Gell- 
Mann (flavor) matrices and A = a/2/31, with 1 being the unit matrix. The fourth term 
is the 't Hooft interaction. [Note that in the notation used in Ref. [2], G s and G v replace 
G s /2 and G v /2 of Eq. (5).] 

In order to make contract with the results of lattice simulations [3-6] we use the model 
with the number of flavors, Nf =1. Therefore, the A* matrices in Eq. (5) may be replaced 
by unity. We then have used 



C = m-rn )q+^[(qq) 2 + (q il5 q) 2 ] 

- ^[{qi^Y + (9757^) 2 ] 



(6) 



in order to calculate the hadronic current correlation functions in earlier work [7-9]. Spectral 
functions obtained using lattice QCD and MEM are shown in Fig. 1 and 2 [3] . 

We note that the unpolarized total cross section for a resonance of mass rriR in the 
s-channel may be written 



—E 
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(7) 



Svr (P 2 -m^ + m|r|' 

where s = 4E 2 in the center-of-mass. The more general result is obtained by using the 
expression 



G, 



9 



so that Eq. (7) becomes 



1 - G s J(P 2 ) 



N 9 
a tot = —E 2 

671 



P 2 -m 2 R + im R TR 



(8) 



G, 



1 - G S J(P 2 



(9) 
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FIG. 1: The spectral functions cj/lo 2 for pseudoscalar states obtained by MEM are shown [3]. The 
solid line is for T/T c = 1.5 and the dashed line is for T/T c = 3.0. The second peak is lattice artifact 
[3]. 
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FIG. 2: The spectral functions a/u 2 for vector states obtained by MEM are shown [3]. See the 
caption of Fig. 1. The second peak is lattice artifact [3]. 
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FIG. 3: The dotted line shown the result for T = 350 MeV obtained in Ref. [2]. (See Fig. 3 of 
Ref. [2] for this result and for results at other temperatures.) The solid curve presents the results 
calculated using our model at T = 350 MeV [7-9] . 




FIG. 4: Quark mean-free path, \{p2), shown for different densities of antiquaries as a function of 
the quark momentum p2- (See Ref. [10] and Eq. (8) of the present work for the parametrization of 
the antiquark densities used in our calculation.) Here, \i = 1.1 GeV (dotted curve), ji = 1.3 GeV 
(dashed curve) and \i = 1.5 GeV (solid curve). 
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where N is a statistical factor which we take to be N = ^2j(2J + 1) = 8, since we consider 
a sum over scalar, pseudoscalar, vector and axial-vector resonances in our model. The 
calculation of J{P 2 ) is discussed in great detail in our earlier work [7-9] and we do not 
repeat that discussion here. (In that work we have used a Gaussian regulator rather than 
the sharp cutoff usually used in the case of the NJL model.) 

The result for the cross section u + u — > u + u obtained in Ref. [2] is shown as a dotted 
line in Fig. 3. The solid line shows the result obtained in our model [7-9]. 

The application of our model to calculate the quark mean-free path is described in Ref. [10] 
and some results are shown in Fig. 4 for different densities of antiquarks in the plasma [10]. 
That distribution is parametrized as in Eq. (2.7) of Ref. [10] with 

= expWH) -„] + !• (10) 

(3 = 1/T and E{p[) = [pi 2 + m 2 } 1 / 2 . (See the caption of Fig. 4.) 

The main difference between our model and that of Ref. [2] is that we recognize that the 
system is deconfmed above the critical temperature. The resonances calculated in Ref. [2] 
using the NJL model that play the major role in the results of that work are not present in 
the deconfmed plasma. In our work we have used the resonances found in the deconfmed 
phase when using lattice QCD and the maximum entropy method. Our work is closer in 
spirit to that of the Stony Brook group [11-13] who also make reference to the QCD lattice 
data in their discussion of the quite small quark mean-free paths. 
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